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Location of the Multicritical Point for the Ising Spin Glasses on the 
Triangular and Hexagonal Lattices 
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A conjecture is given for the exact location of the multicritical point in the phase diagram 
of the ± J Ising model on the triangular lattice. The result p c = 0.8358058 agrees well with a 
recent numerical estimate. From this value, it is possible to derive a comparable conjecture 
for the exact location of the multicritical point for the hexagonal lattice, p c — 0.9327041, 
again in excellent agreement with a numerical study. The method is a variant of duality 
transformation to relate the triangular lattice directly with its dual triangular lattice without 
recourse to the hexagonal lattice, in conjunction with the replica method. 
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1. Introduction 

Properties of finite-dimensional spin glasses are still under active current investigations 
after thirty years since the proposal of basic models. 1,2 ) Among issues are the existence and 
critical properties of spin glass transition, low-temperature slow dynamics, and competition 
between the spin glass and conventional phases. 3 ) 

Determination of the structure of the phase diagram belongs to this last class of problems. 
Recent developments of analytical theory for this purpose, 4 ~ 7 ** namely a conjecture on the 
exact location of the multicritical point, have opened a new perspective. The exact value of 
the multicritical point, in addition to its intrinsic interest as one of the rare exact results for 
finite-dimensional spin glasses, greatly facilitates precise determination of critical exponents 
around the multicritical point in numerical studies. 

The theory to derive the conjectured exact location of the multicritical point has made use 
of the replica method in conjunction with duality transformation. The latter aspect restricts 
the direct application of the method to self-dual lattices. It has not been possible to predict 
the location of the multicritical point for systems on the triangular and hexagonal lattices 
although a relation between these mutually-dual cases has been given. 7 ) 

In the present paper we use a variant of duality transformation to derive a conjecture for 
the ± J Ising model on the triangular lattice. The present type of duality allows us to directly 
map the triangular lattice to another triangular lattice without recourse to the hexagonal 
lattice or the star-triangle transformation adopted in the conventional approach. This aspect 
is particularly useful to treat the present disordered system under the replica formalism as 
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Fig. 1. Triangular lattice and its dual (shown dashed). The up-pointing triangle surrounded by the 
variables fa, fa, fa is transformed into the down-pointing triangle surrounded by fci, &2, £3. 



will be shown below. The result agrees impressively with a recent numerical estimate of high 
precision. This lends additional support to the reliability of our theoretical framework 4-7 ) 
to derive a conjecture on the exact location of the multicritical point for models of finite- 
dimensional spin glasses. 

2. Ferromagnetic system on the triangular lattice 

It will be useful to first review the duality transformation for the non-random Z q model 
on the triangular lattice formulated without explicit recourse to the hexagonal lattice or the 
star-triangle transformation. 8 ) Let us consider the Z q model with an edge Boltzmann factor 
x[fa — fa] for neighbouring sites i and j. The spin variables fa and fa take values from to 
q — 1 (mod q). The function x[-] itself is also defined with mod q. An example is the clock 
model with coupling K, 

x[fa - fa] = exp j if cos (j^-ifa ~ <Aj)) j • (!) 

The Ising model corresponds to the case q = 2. 
The partition function may be written as 

Z = VY1 n X ^ 12 ] X ^ 2 3] X ^3l]^12+^23 + ^3l)n <5 (E^)- ( 2 ) 

{4>ij} a v 

Here the product over A runs over up-pointing triangles shown shaded in Fig. 1 and that 
for y is over unshaded down-pointing triangles. The variable of summation is not written as 
the original fa but in terms of the directed difference faj = fa — fa defined on each bond. 
This is possible if we introduce restrictions represented by the Kronecker deltas (which are 
defined with mod q) as in eq. (2) allocated to all up-pointing and down-pointing triangles. For 
instance, <^i2(= fa — fa)-, faz{= fa — fa), and fai(= fa — fa) are not independent but satisfy 
^12 + 023 + 031 = (mod q), where 1,2, and 3 are sites around the unit triangle as indicated 
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in Fig. 1. The overall factor q on the right hand side of eq. (2) reflects the invariance of the 
system under the uniform change <pi — ► <pi + I (Vi, < / < q — 1). 

It is convenient to Fourier-transform the Kronecker deltas for down-pointing triangles 
and allocate the resulting exponential factors to the edges of three neighbouring up-pointing 
triangles. Then the partition function can be written only in terms of a product over up- 
pointing triangles: 



1 f 2m 

-^[012,023, 03l]<5(012 + 023 + 03l)exp | — (fci023 + /c 2 031 + h(f) 12 ] 



{ki}{4>ij} A 

(3) 

where A[<f>i 2 , 023, hi] = z[0i2M023M03i]- 

Now let us regard the product over up-pointing triangles in eq. (3) as a product over 
down-pointing triangles overlaying the original up-pointing triangles as shown dashed in Fig. 
1. This viewpoint allows us to regard the quantity in the square brackets of eq. (3) as the 
Boltzmann factor for the unit triangle (to be called the face Boltzmann factor hereafter) of 
the dual triangular lattice composed of overlaying down-pointing triangles: 

{hi} V* 

where 

1 q ~ l 

A* [k u , k 23 ,k 31 ] = - Y) 4012, 023, 03l]<5(012 + 023 + 031 ) 
q 

4>12,4>23,4>31=0 

( 2iri 1 
X exp j— — (fci023 + &2031 + h(f) 12 )> . (5) 

Here we have used the fact that the right hand side is a function of the differences ki — kj = 
kij ({ij) = (12), (23), (31)) due to the constraint 4>i 2 + ^23 + ^31 = 0. 

This is a duality relation which exchanges the original model on the triangular lattice with 
a dual system on the dual triangular lattice. A*[ki 2 ,k 2 3,k 3 i] represents the face Boltzmann 
factor of the dual system, which is the function of the differences between the nearest neighbor 
sites on the unit triangles similarly to face Boltzmann factor ^4[0i2, 023, 03i] of the original 
system. 

It is easy to verify that the usual duality relation for the triangular lattice emerges from the 
present formulation. As an example, the ferromagnetic Ising model on the triangular lattice 
has the following face Boltzmann factors: 

40,0,0] =e 3K , A[l, 1,0] =41,0,1] =A[0, 1,1] =e' K ,--- , (6) 

where A[0, 0, 0] is for the all-parallel neighbouring spin configuration for three edges of a unit 
triangle, and A [1,1,0] is for two antiparallel pairs and a single parallel pair around a unit 
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triangle. The dual are, according to eq. (5), 

A* [0, 0, 0] = \ {A[0, 0, 0] + 3A[1, 1,0]}, A* [1, 1, 0] = ^ {A[0, 0, 0] - A[l, 1,0]}. (7) 
It then follows that 

-4K* = A*[1,1,0] 1-e- 4 ^ 

~ A* [0,0,0] l + 3e- 4 ^' { ) 

This formula is equivalent to the expression obtained by the ordinary duality, which relates 
the triangular lattice to the hexagonal lattice, followed by the star-triangle transformation: 

(l + 3e- 4 ^)(l + 3e- 4 ^*) =4. (9) 

3. Replicated system 

It is straightforward to generalize the formulation of the previous section to the spin glass 
model using the replica method. The duality relation for the face Boltzmann factor of the 
replicated system is 



A*m 2 },{k% 3 },{k%i}} = -JrE 



Q 



J] sift + # + C) 



.0=1 



x ex P |^ X>?C + + j , (io) 

where a is the replica index running from 1 to n, {kfj} denotes the set {kjj,kfj,--- ,kfj}, 
and similarly for {</>"} etc. The variables 4>%, 0" correspond to 4>u, 023, 031 in eq. (3). The 
original face Boltzmann factor is the product of three edge Boltzmann factors 

A [{<Pa}, {<Pb}, {€}] = X4>l-n ■ ■ ( H ) 

where x^ 1 ■■■</>« is the averaged edge Boltzmann factor 

9-1 

z=o 

Here pi is the probability that the relative value of neighbouring spin variables is shifted by I. 
A simple example is the ±J Ising model (q = 2), in which po = p (ferromagnetic interaction) 
and pi = 1 — p (antiferromagnetic interaction). 

The average of the replicated partition function Z n is a function of face Boltzmann factors 
for various values of 0's. The triangular-triangular duality relation is then written as 5 ~ 7 ) 

Z n (A[{0}, {0}, {0}], • • • ) = cZ n (A*[{0}, {0}, {0}], • • • ), (13) 

where c is a trivial constant and {0} denotes the set of n 0's. Since eq. (13) is a duality 
relation for a multivariable function, it is in general impossible to identify the singularity of 
the system with the fixed point of the duality transformation. Nevertheless, it has been firmly 
established in simpler cases (such as the square lattice) that the location of the multicritical 
point in the phase diagram of spin glasses can be predicted very accurately, possibly exactly, by 
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using the fixed-point condition of the principal Boltzmann factor for all-parallel configuration 
jO}. 4-7 -* We therefore try the ansatz also for the triangular lattice that the exact location of 
the multicritical point of the replicated system is given by the fixed-point condition of the 
principal face Boltzmann factor: 



A[{0},{0},{0}]=A* [{0},{0},{0}], 



(14) 



combined with the Nishimori line (NL) condition, on which the multicritical point is expected 
to lie. 9 ' 10 ) 

For simplicity, we restrict ourselves to the ± J Ising model hereafter. Then the NL condition 
(1 — p)/p, where p is the probability of ferromagnetic interaction. The original face 



is e 



-IK _ 



Boltzmann factor is a simple product of three edge Boltzmann factors, 

A[{0},{0},{0}]= X 3 

with 

e (n+l)K + e -(n+l)K 



Xo = pe nK + (1 - P)e 



-nK 



e K + e~ K 



(15) 



(16) 



The dual Boltzmann factor A*[{0}, {0}, {0}] needs a more elaborate treatment. The con- 
straint of Kronecker delta in eq. (10) may be expressed as 



^^(C+C + C) 



a=l 



(17) 



n«+c+o=2- n 2 ex p 

a=l {r a =0,l} 

The face Boltzmann factor A [{</>„}, {4>c}] * s the product of three edge Boltzmann fac- 
tors, each of which may be written as, on the NL, 7 ) 



p exp 
1 



+ (1 - p) exp 



Q=l 



2 cosh K 



E exp 



J?o=±l 



a=l 



(18) 



Using eqs. (17) and (18), eq. (10) can be rewritten as, for the principal Boltzmann factor with 
all kf = 0, 



A* [{0},{0},{0}] 



4™ (2 cosh K) 3 



EEE ex P 



r\ t 



7TZJ> Q (C + C + C) 



+ K( Va + i lb + T] C ) + K Va 2(1 - 2<ft) + K Vb ^(l " + K Vc 2(1 " 2#) 



(19) 



The right hand side of this equation can be evaluated explicitly as shown in the Appendix. 
The result is 



A* [{0},{0},{0}] = 4- n (2coshK) 3n - 3 

x [(e 3 ^ + 3e~ K )(l + tanh 3 K) n + (3e^ + e" 3 ^)(l - tanh 3 K) 1 



(20) 
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The prescription (14) for the multicritical point is therefore 

cosh 3 (n+ 1)K _ (2 cosh K) 3n 
cosh 3 K ~ 4 n (2coshK) 3 

x [(e 3 ^ + 3e- K )(l + tanh 3 K) n + (3e K + e~ 3K )(l - tanh 3 K) n ] . (21) 
4. Multicritical point 

The conjecture (21) for the exact location of the multicritical point can be verified for 
n = 1, 2 and oo since these cases can be treated directly without using the above formulation. 

The case n = 1 is an annealed system and the problem can be solved explicitly. It is easy 
to show that the annealed ± J Ising model is equivalent to the ferromagnetic Ising model with 
effective coupling K satisfying 

tanh K = (2p - 1) tanh K. (22) 

If we insert the transition point of the ferromagnetic Ising model on the triangular lattice 
e 4R = 3, eq. (22) reads 

(2p- 1) tanh if = 2- V3. (23) 

This formula represents the exact phase boundary for the annealed system. Under the NL 
condition, it is straightforward to verify that this expression agrees with the conjectured 
multicritical point of eq. (21). It is indeed possible to show further that the whole phase 
boundary of eq. (23) can be derived by evaluating A* [0, 0, 0] directly for n = 1 for arbitrary p 
and K, giving 2,4* [0, 0, 0] = Xo + %Xi Xo with X o = pe K + (1 - p)e~ K , Xi = P?~ K + (1 - P> K , 
and using the condition A*[0, 0, 0] = A[0, 0, 0]. 

When n = 2, a direct evaluation of the edge Boltzmann factor reveals that the system on 
the NL is a four-state Potts model with effective coupling K satisfying 5, 7 ) 

e * = e 2K -l + e~ 2K . (24) 

Since the transition point of the non-random four-state Potts model on the triangular lattice 
is given by e K = 2, 8 ) the (multi) critical point of the n = 2 system is specified by the relation 

e 2K _ ! + e -2K = 2 (25) 

Equation (21) with n = 2 also gives this same expression, which confirms validity of our 
conjecture in the present case as well. 

The limit n — ► oo can be analyzed as follows. 5 ) The average of the replicated partition 
function 

[Z n U = [e-"^] av , (26) 

where [■ ■ ■ ] av denotes the configurational average, is dominated in the limit n — > oo by contri- 
butions from bond configurations with the smallest value of the free energy F. It is expected 
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that the bond configurations without frustration (i.e., ferromagnetic Ising model and its gauge 
equivalents) have the smallest free energy, and therefore we may reasonably expect that the 
n — > oo systems is described by the non-random model. Thus the critical point is given by 
e 4K = 3. It is straightforward to check that eq. (21) reduces to the same equation in the limit 
n — ► oo. 

These analyses give us a good motivation to apply eq. (21) to the quenched limit n — > 0. 
Expanding eq. (21) around n = 0, we find, from the coefficients of terms linear in n, 

3Kt<mhK = 3 log (2 cosh K) - log 4 

pZK , o „-K o p K , -3K 

+ (^W los(1 + tanh K) + (2^W los(1 - taDh K) - (2?) 

or, in terms of p, using e~ 2K = (1 —p)/p, 

2p 2 (3 - 2p) logp + 2(1 - pf{l + 2p) log(l -p)+ log 2 

= p(4p 2 - 6p + 3) log(4p 2 - 6p + 3) + (1 - p)(4p 2 - 2p + 1) log(4p 2 - 2p + 1). (28) 

Equation (28) is our conjecture for the exact location of the multicritical point p c of the ± J 
Ising model on the triangular lattice. This gives p c = 0.8358058, which agrees well with a 
recent high-precision numerical estimate, 0. 8355(5). n ) 

If we further use the conjecture 7 ^ H(p c ) + H(p' c ) = 1, where H(p) is the binary entropy 
— plog 2 p — (1 — p) log2(l — p), to relate this p c with that for the hexagonal lattice p' c , we find 
p' c = 0.9327041. Again, the numerical result 0.9325(5) n ) is very close to this conclusion. 

5. Conclusion 

To summarize, we have formulated the duality transformation of the replicated random 
system on the triangular lattice, which brings the triangular lattice to a dual triangular lattice 
without recourse to the hexagonal lattice. The result was used to predict the exact location 
of the multicritical point of the ± J Ising model on the triangular lattice. Correctness of our 
theory has been confirmed in directly solvable cases of n = 1,2 and oo. Application to the 
quenched limit n — > yielded a value in impressive agreement with a numerical estimate. 

The status of our result for the quenched system, eq. (28), is a conjecture for the exact 
solution. It is difficult at present to prove this formula rigorously. This is the same situation as 
in cases for other lattices and models. 4 ~ 7 ) We nevertheless expect that such a proof should be 
eventually possible since a single unified theoretical framework always gives results in excellent 
agreement with independent numerical estimations for a wide range of systems. Further efforts 
toward a formal proof are required. 
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Appendix 

In this Appendix we evaluate eq. (19) to give eq. (20). Let us denote 
4 n (2cosh A) 3 A*[{0}, {0}, {0}] as A. The sums over a in the exponent of eq. (19) can be 
expressed as the product over a: 



V 

1 1 1 



n 



X 

a=l 



j a =o \<t>%=o <t>%=o 0«=o 

By performing the sums over 4> and r for each replica, we find 



.(29) 



A = e R ( n+1 ^+v b +Vc) Yl [(1 + e" 2 ^)(l + e~ 2 ^)(l + e~ 2K ^) 

Va,Vb>Vc = ±l a=l 

+ (1 _ e ~ 2K Va ) ( l _ e ~2K Vb ^ x _ e -2K Vc _ ^ 

It is straightforward to write down the eight terms appearing in the above sum over r] a ,Vb, ijc 
to yield 

~ A = e 3K(n+l) [ (1+e -2K ) 3 + (1 _ e -2X ) 3]- 

+ 3e K ^ [(1 + e" 2 ^) 2 (l + e 2 *) + (1 - e" 2 ^) 2 (l - e 2 *)f 
+ 3e- K ^ [(1 + e~ 2 ^)(l + e 2 ^) 2 + (1 - e - 2 *)(l - e 2 ^) 2 f 

+ e -3AT(n+l) [ (1 + e 2K )3 + (1 _ ^3]" ^ ^ 

which is further simplified into 

A = e 3K [(2 cosh A) 3 + (2 sinh A) 3 ] n 
+ 3e K [(2 cosh A) 3 - (2 sinh A) 3 ] n 
+ 3e~ K [(2 cosh A) 3 + (2 sinh A) 3 ] n 
+ e" 3 ^ [(2 cosh A) 3 - (2 sinh A) 3 ] n 
= (2coshA) 3n 

x [(e 3K + 3e~ K )(l + tanh 3 K) n + (3e K + e~ 3K )(l - tanh 3 K) n ] . (32) 

This is eq. (20). 
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